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Abstract. - We investigate numerically and on the scaling level the adsorption of a charged 
colloid on an oppositely charged flexible membrane. We show that the long ranged character 
of the electrostatic interaction leads to a wrapping reentrance of the complex as the salt con- 
centration is varied. The membrane wrapping depends on the size of the colloid and on the 
salt concentration and only for intermediate salt concentration and colloid sizes we find full 
wrapping. From the scaling model we derive simple relations for the phase boundaries between 
the different states of the complex, which agree well with the numerical minimization of the 
free energy. 



When a charged spherical particle adsorbs on an oppositely charged flexible membrane, 
the membrane deforms and wraps around the particle. This mechanism plays an important 
role in many physical and biological processes. Introduction of genes into cells via cationic 
lipid-DNA complexes has been one important theme in biological research over the last few 
decades. Recent investigations show that endocytosis of lipid-DNA-complexes is triggered 
by electrostatic interactions with the cell membrane 1,2,3 . Moreover, some viruses rely on 
internalization via endocytic uptake, where non-specific electrostatic interactions (and not a 
particular structural motif) seem to play a crucial role 21[SJ|niE|. In all cases the transfection 
efficiency depends on the size of the interacting particle 8 9 , on the charge densities of the cell 
and particle surfaces and on the salt concentration [2, 3, 4, 8, l(]j. But also in many physical 
experiments wrapping of particles by membranes is important m^OEII- Motivated by 
these examples, we investigate in this paper the role of the electrostatic interaction between a 
charged colloidal particle and an oppositely charged membrane in the adsorption and wrapping 
process. The adsorption of a neutral spherical particle on a neutral fluid membrane due to a 
local adhesion potential has been studied theoretically already in discussing the role of 
the elastic parameters of the membrane in the wrapping process. Furthermore, the wrapping 
of a colloid by a membrane is related to the budding of vesicles, which was subject to many 
studies fEHEj. In the electrostatic case considered in this paper, the attraction between 
the colloid and the membrane as well as the repulsion between the membrane segments are 
shown to play a crucial role. We find that it is the long-ranged character of the electrostatic 
interaction, which leads to a reentrance phenomena depending on the salt concentration: In 
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case of low salt concentration or small colloids, the electrostatic interaction is barely screened 
and due to membrane-membrane repulsion the membrane only slightly wraps around the 
colloid. By increasing the salt concentration or the size of the colloid the screening of the 
membrane- membrane repulsion leads to a wrapping of the colloid by the membrane. Increasing 
the salt concentration or the size further, the electrostatic attraction between the colloid 
and the membrane becomes too weak compared with the mechanical bending energy of the 
membrane and again the membrane is only slightly wrapped. Therefore, the membrane wraps 
the colloid optimally only for intermediate salt concentrations or colloid sizes. This is similar to 
the complexation behavior of a polyelectrolyte and an oppositely charged sphere as discussed 
inCZj. 

In our model, the membrane is described as a two dimensional elastic surface with homoge- 
neous charge density <r, bending rigidity Kq and surface tension 7; the colloid is approximated 
by a homogeneously charged sphere with radius R and charge Z. We describe the electro- 
static interactions by Debye-Hiickel (DH) potentials. The simplest free energy to describe the 
system considered here, in units of /c#T, is given by the sum of the following three parts [T£] : 

F mech = ^ J d 2 p^(p)(AX(p)f+ leff J d 2 p^(p) (1) 



l B Zoe" R / - / _exp{-/t|X(p)|} 



J d 2 p^J) 



1 + kR J r^v, |X(p)| 



(3) 



The vector X(p) describes the position of the membrane in space, where p = (p\, P2) are 
the internal coordinates of the membrane. g(p) is the determinant of the induced metric field 
defined by gij = <9iX<9jX and A = g~?di (y/gg^dj} is the Laplace-Beltrami operator on the 
surface. The purely mechanical (elastic) contribution to the free energy F mec h consists of the 
bending energy (first term) and the surface tension (second term). The repulsive electrostatic 
interaction between the membrane segments is given by F rep and the electrostatic membrane- 
colloid attraction is denoted by F a ur- The inverse Debye-Hiickel screening length is given 
by k := ^/8tt£bc (c is the salt concentration), while the Bjerrum length £3 •= e 2 /ATreosksT 
measures the distance at which two elementary charges interact with thermal energy ksT. 
The membrane is fixed at a circular frame of radius L connected to an area reservoir. In this 
way the projected area of the membrane remains constant as the membrane wraps around 
the colloid. In particular, this means that also the total charge of the system is not preserved, 
which generates an electrostatic tension j e i = —ttIbct 2 /^- Therefore the effective tension 
jeff = 7 + 7d appearing in eq.JIJ consists of the sum of the bare tension 7 and the electrostatic 
tension 7^. 

We analyze the free energy of the system, F = F mec h + F rep + F attr , using two approaches. 
First, we minimize the free energy numerically and, second, we investigate the free energy in 
a simplified geometry on a scaling level. We distinguish three different states of the complex, 
which can be defined as the membrane wraps around the colloid using the wrapping height 
h defined by the point at which the membrane detaches from the colloid (see fig^c)): i) 
the point contact state: the membrane touches the sphere only in one point (h = 0), ii) the 
touching state: membrane touches the colloid over a finite area, but less then half of the 
colloid (R > h > 0) iii) the wrapped state: the membrane covers more then half of the colloid 
surface (2R > h > R). Defining the wrapping threshold at h = R and not at h = 2R (as 
was done e.g. in 14 ) is necessary in order to avoid the singularity at h = 2R caused by the 
touching of charged surfaces in the numerical minimization. In what follows, we rescale the 
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free energy F : = F jkBT2i\K^ and use rescaled variables: Z := Z/AiraR 2 , a := o IbR? I 
7 := syR 2 /Ko and kR. 

In the numerical minimization of the free energy we use a frame with radius L = 2.5ttR 
and model the hard-core repulsion between membrane and colloid by a shifted and truncated 
Lennard- Jones Potential, v LJ (r) = g[(R/\r\) 12 - 2(R/\r\) 6 + 1] for |r| < R and v LJ (R) = 
otherwise. We choose the parameter g in such a way that the equilibrium membrane-sphere 
separation is 0.999.R. Beside the constraint of constant projected area we impose, if not stated 
otherwise, the boundary condition that the membrane becomes planar at the edge. Fig[IJa) 
shows for three different colloid charges Z the numerically obtained membrane shapes for 
kR = 50 (high salt/large colloid limit) and vanishing surface tension 7. For Z = 7.82 the 
membrane is in the touching state. By elevating Z the membrane bends further around the 
colloid and is for Z = 8.38 right at the transition point between the touching and the wrapped 
state, while for Z = 8.54 the membrane is in the wrapped state. Due to the strong screening 
of the electrostatic interaction the dominating energy contribution stems from the mechanical 
bending energy and the membrane wraps into a catenoid-like shape (surface of revolution 
of the curve x(s) « ln(s) for arc- length s ^> 1 [El)? reducing in this way its elastic bending 
energy. In fig^b) we show the unwrapping of the membrane by increasing the surface tension 
for kR = 10 and a = 1.0 (the numbers denote 7). For zero surface tension the membrane 
is in the wrapped state. Increasing the surface tension the membrane unwraps and is for 
7 = 1 right at the transition point to the touching state. By increasing the surface tension 
further to 7 = 2.0 the membrane goes into the touching state. In case of small nR (low 
salt concentration/small colloids), the electrostatic interactions are dominant and thus govern 
the shape of the membrane. In this case the long-ranged electrostatic attraction pulls the 
membrane closer to the colloid, producing the volcano-like shape presented in ngOfd) for 
kR = 0.1 and a = 1.0 (the numbers denote Z). Comparison of the wrapped states in figEJa) 
and figEId) demonstrates clearly the different wrapping structure of the membrane for small 
and large kR, respectively. The shapes in figGfc) are obtained without imposing a boundary 
condition for the membrane at the edge in order to compare the wrapping shapes (symbols) 
of the membrane for different nR and Z with a catenoid (solid line) touching the colloid at 
h = R. The rescaled colloid charge Z was adjusted in such a way that the wrapping height of 
the membrane is in all cases the same (h = R). It can be clearly seen that the deviation from 
the catenoidal-shape becomes pronounced as nR decreases, since the range of the electrostatic 
interaction increases. The calculation of membrane shapes for different parameters allows us 
to measure the complexation diagrams presented in figEl Before we discuss these diagrams we 
calculate the phase boundaries between the different states of the complex using a simplified 
scaling model. 

In the geometry shown in ng^f e) the calculation of the free energy is analytically tractable, 
yielding the free energy in form of a Landau function with the rescaled wrapping height 
h := h/R playing the role of the order parameter, from which we can extract simple scaling 
relations for the phase boundaries. The mechanical energy and the electrostatic membrane- 
colloid attraction read: 




(5) 



(4) 



For the electrostatic repulsion we find in the limit kR <C 1: 
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Figure 1 - Figure a, b, c and d show wrapping shapes of the membrane, a): High salt/large colloid 
limit (kR — 50) for 7 = and a = 1.0. b): Unwrapping of the membrane by increasing the surface 
tension is shown (cr = 1.0, kR = 10.0 and Z = 2.90). c): Increasing deviation from the cateniodal 
shape of the membrane shapes by decreasing kR. The numerically obtained membrane surfaces 
(symbols) as well as the catenoid (solid line) detach from the sphere all at the same wrapping height 
(h — R). The parameters are: kR = 45 , 20 , 10 , 2 and Z = 7.46 , 2.8 , 1.1 , 0.52 (squares, from bottom 
to top), k,R = 0.5, 0.1 and Z — 1.34, 2.22 (triangles, from top to bottom), a = 1 and 7 = in all 
cases, d): Low salt/small colloid limit (kR = 0.1) for 7 = and a = 1.0. e): For the scaling model 
we subdivide the membrane into a spherical segment, which touches the sphere and a planar ring. 



and in the limit kR ^> 1: 



TP — 



kR 2 



l§^~ h 3/2 + ^ _ n ~2 _ 27T KRa 2 ^ h 2 + G(kR- 2 , h^ 2 ) : Vh < V2/kR 



4V2kR 2 



~ h 3/2 



1QV2kR 2 



h h ' 2 + 0(kR- 3 , h 3 ):Vh> V2/kR 



(7) 



Using the equations Q-Q) we can analyze the transitions between the different membrane 
states in the following way. The touching transition occurs when the coefficient in front of h 
vanishes, yielding the exact result for the touching threshold: 

" KR+1 (8) 



Z t = 



2wa 2 



In the high salt limit the free energy given by the sum of the equations Q, (0) and Q 
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Figure 2 - Complexation phase diagrams for the membrane-colloid system as a function of Z and 
kR on a logarithmic scale. In fig. (a) the surface tension is zero. The lower line denotes the touching 
transition, eq.©, and the upper line the wrapping transition, eq.@ and ea. (fTUJ) . The difference 
between the touching and the wrapping threshold is ~ {kK)~ x according to our scaling theory. Thus, 
the touching phase vanishes for large kR. Fig.(b) shows the touching and the wrapping transition 
for finite surface tension, 7 = 10. In both figures the stars denote the results for touching threshold 
and the circles the results for wrapping threshold obtained from the full model. The rescaled charge 
density on the membrane is in all cases given by 5 — 1.0. 



exhibits a secondary minimum rendering the wrapping transition discontinuous. Note that 
this discontinuity is of electrostatic origin and is different from the one found in ^3]. We 
determine the wrapping threshold Z w by resolving the system of equations F(h\) = FQ12) 
and dj l F{h2) = for Z subject to the constraint Y12 < 2, where hi is the location of first and 
h2 is the location of the second minimum. In the limit kR <C 1 the threshold reads: 



1 2 + 7 2y"2 1 



AkR 47TO- 2 7T 4 

and in case of kR > 1 we find: 



(D(kR) : . (9) 



The resulting complexation diagram in the Z j kR plane is presented in fig El where the lines 
denote the scaling results eq. (|5|)-(|10[) and the symbols refer to the numerical solution of the 
full model. For kR small the touching transition depends in leading order only on the rescaled 
charge density <r, while in the high salt limit the touching threshold is Z t ~ nR/d 2 . If the 
surface tension is zero (figEfa)), the touching threshold differs from the wrapping threshold for 
kR ^> 1 by a term of O ((fti?) -1 ) . It follows from eq.(|EJ) and ec fTHj) that in the limit kR — > 00 
the membrane goes for ZoIbi 1 k = 2Ko directly from the point contact state with h = to the 
wrapping state with h > 1, skipping the touching region 20,21 j. This limit corresponds to 
the case of a local adhesion energy between a neutral colloid and neutral membrane and was 
discussed in [T3) . We find very good agreement between the phase boundaries derived from 
the scaling model and the values for touching and the wrapping threshold obtained from the 
full model, demonstrating that the simplified model captures the main physical properties of 
the system. 

Starting from the touching state in fig El an d fixing Z the complex goes by increasing the 
salt concentration into the wrapping state and by increasing the salt concentration further 
finally leaves the wrapping state and goes into the touching state again. This reentrance 
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a) b) 




Figure 3 - a): The wrapping height h as a function of kR. The line denote h as obtained from 
the scaling model (eq.0)-eq.Q), while the solid circles denote the corresponding results obtained by 
numerical minimization of the full free energy. The parameters are: a — 1.0, Z — 2.2 and 7 = 5.0. 
b): the complexation phase diagram for the membrane-colloid system as a function of Z and a is 
shown on a logarithmic scale (kR = 10, 7 = 5). The lower line denotes the touching transition and 
the upper line shows the wrapping transition. The stars symbolize the touching threshold and the 
solid circles the wrapping threshold as obtained from the numerical solution of the full model. 



behavior of the complex can also be seen in figCHa), where we plot the wrapping height h : 
which minimizes the scaling free energy, as a function of kR. The solid line results from the 
scaling free energy, eq.Q-Q with the constraint h < 2 and the solid circles show the results 
obtained from the full model. The physical mechanism for this wrapping/unwrapping is the 
following. If kR is small, the electrostatic interaction is long ranged and the free membrane 
area is strongly self repelling. As the membrane bends around the colloid, charged membrane 
area is pulled closer together and the resulting energy cost has to be paid for by increasing 
the membrane-colloid attraction. In particular, this means that for zero salt concentration 
(kR = 0), i.e. in the case of pure Coulomb interaction, the membrane only slightly bends 
around the colloid, no matter how large Z is. On the other hand, if kR is large the electrostatic 
membrane-colloid attraction is screened and in order to overcome the bending energy of the 
membrane the rescaled colloid charge has to be increased again. 

When the bending rigidity Ko is small and/or the surface charge density on the mem- 
brane is large the parameter a := cr^ £bR 3 / 'Ko becomes large. In this case the electrostatic 
interaction dominates the behavior of the complex. The limit a — > 00 corresponds to the 
case of vanishing bending rigidity and even an infinitesimal small charge on the sphere will 
result in a bending of the membrane, i.e. Z t — > as follows from eq.®. On the opposite, the 
limit d — > corresponds to infinitely large bending rigidity Kq (stiff wall) or vanishing surface 
charge density a on the membrane and the membrane remains planar as long as Za — > 0, 
which also follows from eq.Q. In figCJb) we show the complexation diagram in the Z/a-plane 
(kR = 10, 7 = 5). The symbols denote the results from the full model and the solid lines the 
results from the scaling model eq. (JSJ) and ea.([TUj). Again, we find good agreement between 
the simple scaling model and the much more complex full model. The wrapping threshold 
diverges as Z w ~ a~ 2 for a — > 0, while for a — > 00 we find Z w ~ 1/nR, which follows from 

In this paper we have studied the role of the electrostatic interaction between a charged 
colloid and an oppositely charged fluid membrane in the wrapping process. For nR <C 1 the 
membrane wraps into a volcano-like shape, while in the opposite limit nR ^> 1 the membrane 
wraps into a catenoid-like shape. We find a reentrance behavior in the wrapping process 
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depending on kR. For small as well as for large kR the membrane wraps the colloid only 
for large colloid charges. From the scaling model we derive simple relations for the phase 
boundaries of the complex, which agree well with our results from the more complex full 
model. The electrostatic repulsion between the membrane segments plays a crucial role in 
the adsorption process of charged colloid on an oppositely charged membrane. But also for 
the adhesion of neutral colloids at charged membranes in the case of a non-electrostatic local 
adhesion energy the membrane self repulsion is important. Since most bio- membranes are 
negatively charged the electrostatic mechanism studied in this paper plays certainly a role in 
many biological systems. 
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